High spin proton emitters in odd-odd nuclei and shape changes 
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Abstract 



We present a formalism to describe proton emission from odd-odd nuclei 
based on a scattering like approach. Special emphasis is given to the case 
of transitions between states with different deformations. As an example we 
estimate the proton half-life of the odd-odd nucleus ^®Cu. Our calculations 
show that the change of deformation in the decay process has a significant 
influence on the half-life. In addition, the angular momentum coupling of 
proton and neutron orbitals can result in an important ii'-hindrance of the 
decay. To fully account for the observed half- life in ^^Cu, one has to consider 
a shape mixing in the final state. 



PACS numbers: 23.50. +z, 24.10.Eq, 27.40.+Z 
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The detection and study of nuclei on the proton drip hne are being pursued intensely at 
present Many spherical as well as deformed nuclei in the region 50 < Z < 82 have been 
shown to be proton emitters . Recently, an interesting example of prompt 

proton decay was reported: the proton decay of the medium mass odd-odd nucleus ^^Cu 
which includes a transition from a highly deformed state into a near spherical state, thus 
showing a pronounced shape change during the transition from the mother to the daughter 
nucleus pUj] , 

The study of proton emission are mainly based on a scattering like approach 
TT,I2,r3,Tl,rB[ or on the R-matrix formalism |TB| , pT| , p!B| . A comparison of the two approaches 



was recently performed in Refs. |]T^|T9|. These studies were restricted to the case of decay 
from odd-even to even-even nuclei. Only recently odd-odd nuclei have been considered pO[| . 
However, in all cases it has been assumed that the shapes of the mother and the daughter 
nuclei are the same. 

In this letter we present a formalism based on the scattering like approach to study 
proton emission from deformed odd-odd nuclei including transitions between excited states 
as well as different deformations between the mother and daughter nuclei. In addition, we 
investigate the role played by the angular momentum orientation of the odd particle in the 
intrinsic system. As an example we will consider the case of ^^Cu, which life time has been 
established recently pT| . 

We will present in some detail the case corresponding to proton emission processes con- 
necting two axially symmetric deformed nuclei. 



The odd-odd mother wave function in the laboratory frame is given by the ansatz 
1 



•(TfJi 



/2(1 + Sk^a) 

where Q are the Euler angles and the normalised Wigner function 



V'^,^^X^) ^ ^^-^D'^^^^iQ) (2) 

describes the rotation of the core. The prime coordinates (e. g. r^) refers to the intrinsic 
frame. 

The intrinsic wave function of the odd-odd (proton-neutron) system is described by 

^PKA<,<) = ^KA<)^KM) , (3) 

where the initial intrinsic spin projection is Ki = K^^ + K^. The functions $ are standard 
Nilsson wave functions which can be written in terms of spherical components, i. e. 



$^.(r;) = E^^^^3^^.x.(r:), (4) 



, r = 7r,z/. 



In order to evaluate the proton decay width we express the mother wave function in terms 
of the coordinates of the outgoing proton in the laboratory frame. By performing standard 
manipulations the wave function (|l|) becomes a superposition of terms corresponding to the 
proton wave function coupled to the rotational band built on the odd neutron state with 
spin projection K^, i. e. 
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where 
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By using the orthogonahty relation of these functions 

{4>Jfij\4>j'^.i'j') = {Ji, Ki] j, -K^\Jf, K^)^6j^j'^6ii'6jj 



(7) 



the proton partial decay width to a fixed final state Jt becomes 



Tj^ = hv^{Ji, Ki] j, -K^\Jf, hm^\gijK^{r^)\'^ . 

ij 

To obtain the radial parts of the wave functions, i. e. the functions gijK-^ifn) in Eq- 
has to solve the corresponding Schrodinger equation, i. e. 



one 



(9) 



where H^^ is the Hamiltonian describing the motion of the emitted proton with positive 
energy E^^. We assume that the neutron is a spectator. That is, the Hamiltonian H^^ will 
not affect the neutron coordinates. 

One can solve Eq. (j^) by transforming to the intrinsic frame, as in Eq. ([l|). One then 
integrates over the neutron coordinates and all angular coordinates to obtain the system of 
equations corresponding to the radial proton wave functions, i.e. 
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{yijKjVi\yi'j'Kjgi'j'K^ + {yijKjV2\yi'fKj 



Here V2 is given by the non-spherical part of the spin-orbit interaction [|15|1 . 

Our integration procedure using all possible channels confirms that the initial state in the 
odd-odd nucleus has a "stretched" configuration 7igQ/2 ® [0 with an amplitude 0.98. 
For the initial state one therefore has Jj = 2j = 9. As a result the system of equations ( |T0| ) 
practically contains only that component. That is, the outgoing proton escapes as a wave 
containing all angular momenta but the contribution of the configuration = (4,9/2) is 
clearly dominant. Therefore, in our radial equations we match the external solution with 
only that internal wave function. We impose on the radial component gijK-„{f^n) outgoing 
boundary conditions at large distances, i. e. 



'111 



where Fi{krT^) and Gi{kr.,r) are the regular and irregular Coulomb functions, respectively. 
'H^^^^lkr^,) denotes the numerical solution found using backward integration. The wave 
number is given by 
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k = ^2nE^/h 



(12) 



The constant C in Eq. (|TT]) is obtained by matching the external and internal solutions at 
some point R 

CijK.=Vi,KAR)/\ni^^ikR)\ . (13) 

where VijK„{R) is the proton formation amplitude. 

Since lim^^^oo \i^i{kr.^)\'^ + \Gi{kr.^)\'^ = 1 one can write the final expression for the proton 
decay width (||) as 



{J„Kf,j,-K^\Jj,K,)Vi,KAR)/\nl^^kR)\\- . (14) 
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A similar formula was already presented in Ref. pO . 

We will assume that the mother and daughter even-even cores can be described within 
the BCS approach. As already mentioned the odd neutron acts as a spectator. The proton- 
neutron interaction is included in a phenomenological manner by adjusting the proton and 
neutron pairing gaps to their experimental values. With these approximations the formation 
factor V is computed as the internal component multiplied by the amplitude to find the final 
proton Nilsson state (cj^x^) in the initial quasiparticle state (ajj^^), i.e. 

VukAR) = 9tKliR){BCSf\c,,KAjKjBCS.) , (15) 

where QiJkI is the internal solution of the system (|T0[) and \BCSi(^f)) refers to the initial 
(final) ^^Ni core state. One finally obtains 

Vi.kM) = gtKliRWilUK^KABCSflBCS,) . (16) 

Here m^^(t) is the proton BCS amplitude of the final nucleus. The elements of the matrix U 
connect the quasiparticle operators corresponding to the initial and final deformations (for 
details see Appendix (E.3) of Ref. [0]). The overlap integral of the two cores is given by 

/ = (BCSflBCSi) = \det[U{7T)]det[U{u)p^ , (17) 

where ^'"^P denote the single particle Nilsson wave functions. In Eq. (^) we assume that 
the X-value for both the odd-proton and neutron states are conserved during the decay 
process. 

The proton wave function necessary to evaluate the formation amplitude in Eq. (p 



will be calculated by using a Woods-Saxon mean field with universal parametrisation |2S]. 
We adjust the depth of the potential in order to reproduce the Q-value. 

The decay width (|I^), or the half-life T = Ti In2/T, should be independent of the radius 
i? in a region beyond the nuclear surface. This condition is automatically satisfied in our 
case because the internal and external solutions satisfy the same Schrodinger equation. 

To analyse the i^-dependence and to determine the shape of the mother and daughter 
nuclei, we performed total Routhian surface (TRS) calculations, where quadrupole and hex- 
adecapole deformations are minimized with respect to the energy in the rotating frame of 
reference. Pairing correlations are treated selfconsistently by means of the Lipkin-Nogami 
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method p6| , p7| , |28[] . The shape of the mother (initial) nucleus ^^Cu is calculated to be pro- 



late and axially symmetric. The hexadecapole deformation is calculated to be = 0.05 
while the quadrupole deformation parameter yields a value of /?2 = 0.35. The bandhead 
spin of the mother nucleus, assigned to be J = 9, corresponds to a predominantly K = 1/2 
configuration for both protons and neutrons, that is fully rotationally aligned. 

The initial i^'-value of ^^Cu is given by Ki = ± K,y\ yielding a value of Ki = 1, 0. 
If one assumes that the daughter (final) nucleus ^^Ni has the same prolate deformation as 
^^Cu the overlap is unity and the half-live becomes T ^ 0.8 x 10~^^ s, which underestimates 
the corresponding experimental value, i. e. Tg^p ^ 0.5 x 10~^^ s [3T|, by four orders of 
magnitude. 

Actually the nucleus ^^Ni is expected to be spherical. However, our calculations yield 
the shape of the J/ = 9/2"*" state to be oblate deformed, with a deformation value of ~ 
—0.2. This oblate deformed state is approximately 600 keV lower in energy as compared 
to the spherical minimum. We thus expect the ground state of the Jf = 9/2"*" state to 
be oblate deformed and not spherical. Since the moment of inertia at oblate shapes is 
rather small, non collective excitations can compete or mix into the corresponding collective 
rotations. According to the TRS calculations, the oblate collective minimum disappears 
with increasing angular momenta, whereas the non-collective oblate minimum stays stable 
up to high frequencies. The quadrupole deformation lifts the degeneracy of the Jf = 9/2"*" 
multiplet, forming a distinct state with a i^'-value of 9/2 as ground state. For collective 
rotational states one expects a strongly coupled band involving Ml transitions. Indeed, the 
data on ^'''Ni shows a Ml structure above the proposed Jf = 9/2"*" state, distinctly different 
from the rest of the level scheme |]29|. The first excited state with spin 11/2"*" is also ~1 MeV 
lower in energy than the corresponding excitations of the negative parity spherical states. 
Although the structure resembles that of a strongly coupled band, it does not follow a J( J+1) 
pattern, pointing towards non collective excitations, in agreement with the calculations. 

Assuming that the mother nucleus is prolate and the daughter spherical the overlap 
([T7| ) becomes P ^ 3.2 x 10~^, giving a half-life of T ^ 0.24 x 10"^"^ s, which is still two 
orders of magnitude below the experimental value. If one assumes the oblate shape discussed 
above for the daughter nucleus the half life increases further and one obtains P ~ 10~^ and 
T ^ 0.7 X 10~^^ s. The change is not large as compared with the prolate-spherical transition 
assumed above. 

In order to understand the discrepancy with the experimental halflife T, we analysed in 
detail the different elements entering the expression of the decay width. We thus noticed 
that the ii'-values corresponding to the odd proton and neutron which may couple to the 
initial projection Ki, ranging from K=l/2 to K=9/2, can change drastically the values of the 
Clebsch-Gordan coefficient in Eq. ([l^)- In Fig- 1 we show the half-life computed according 
to the decay width given by Eq. (|^) as a function of the final (neutron) spin projection in 
the intrinsic system Kf = K^. The solid lines correspond to decays where the mother and 
daughter nuclei have the same deformation, i. e. where jSi = jSf = 0.35 and the overlap 
factor entering Eq. (P^ ) is virtually unity. In addition to Ki = 1, (upper and middle line) 
we also considered the value Ki = 2 (lower line). In the same figure we show the effect of 
the change in deformation due to the overlap, Eq. (|1^, for prolate and oblate shapes. The 
dashed lines correspond to the case where the daughter nucleus has an oblate deformation, 
i. e. where Pi = 0.35 and Pf = —0.2. 

The oblate deformation has a profound influence on the wave function of the daughter 
nucleus: Whereas for spherical nuclei all K-vsAues are degenerate and equally probable, at 
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oblate shapes the K.,^ = 9/2 orbit of the (79/2 subshell is at the Fermi surface. This imphes, 
according to Eq. (p!^ , that the transition is forbidden since K^^ = 1/2 and K^, = 9/2 would 
result in a vanishing Clebsch Gordan coefficient. Thus, in addition to the shape hindrance 
we encounter a K-hindrance. 

Therefore the only possibility to account for the proton decay of ^^Cu is to assume that 
the initial and final i^'-values are not pure. This is in fact expected, since the Coriolis 
interaction will induce i^'-mixing. To further investigate the role of this interaction, we 
performed additional particle plus rotor calculation for the initial and final nuclei by using 
the Nilsson basis functions given by Eq. (^, i.e. 

=T.d%''^%'^ . (18) 
The odd neutron is not anymore a spectator and its preformation factor is given by 

Our calculation shows that different eigenstates with the same projection give a small con- 
tribution toVy. 

The effect of the Coriolis mixing is two-fold: 

i) The initial and final Ki, will now acquire components of higher and lower values, as 
can be seen in Table 1. From the same Table it becomes clear that indeed the predominant 
component in the initial prolate nucleus has Ky = 1/2. 

ii) The value of corresponding to the band-head will change accordingly allowing 
additional coupling possibilities. The new term arising from the Coriolis mixed states implies 
that the K.„ = ±1/2 can couple to Ki, = =f5/2 yielding a. Ki = 2, 3 state. According to 
Table 1, this is the only common component in both initial and final nuclei. Higher i^-values 
can be possible although with even smaller amplitudes. 

Then, assuming that the daughter nucleus is deformed with jSf = —0.2 the calculated 
half-life for the case of i^'j = 2 and Kf = = 5/2 (see Fig. 1) becomes T ^ 1.8 x 10"^^ s. 
The neutron preformation factor is ~ 10~^. As a result the lifetime acquires the value 
T 1.8 X 10~^ s, which overestimates the experimental value by three orders of magnitude 

Therefore, the experimental value would correspond to a decay in the daughter nucleus 
where the 9/2"*" state has components in the wave function corresponding to oblate and spher- 
ical shapes, respectively. In other words, shape fiuctuations due to the relative /3-softness 
of the final nucleus provides a suitable scenario to explain the actual lifetime. Assuming 
that the spherical and oblate states in ^^Ni are mixed, one obtains a lifetime between the 
spherical T ^ 10"^^ s and oblate prediction T ^ 10~^ s. A simple estimate shows that 
the amplitudes of spherical and oblate components should correspond to Xsph ~ 0.1 and 
^obi ~ 0.9, respectivelly, in order to reproduce the experimental half-life. 

In conclusion we have presented in this letter a formalism to analyse proton decay in 
situations where the initial and final nuclei are in excited states with different deformations. 
We applied the formalism to study the decay of the odd-odd nucleus ^^Cu and found that 
due to the difference of deformations between the mother and daughter nuclei the calculated 
half-life increases by two orders of magnitude, but still remains two orders of magnitude 
below the experimental value. 

Due to the prolate shape of the mother and oblate shape of the daughter nucleus, the 
decay is K-hindered and can proceed only via a very weak component originating from the 
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Coriolis mixing of initial and final states. The half-life in this case becomes three orders 
of magnitude larger than the experimental value. A possible explanation of the measured 
half life is to assume mixing of spherical and oblate shapes in the daughter nucleus, i. e. in 
^^Ni. The spherical component has a small amplitude pa 0.1 but it is essential to describe 
the experimental half-life. 

Given the large changes in life times, associated with the different physical processes 
involved in the decay, implies of course a substantial uncertainty in our analysis. Still, these 
unexpected features indicate that proton decay is a powerful tool to investigate fine details 
of the intrinsic structure of nuclei at the drip line. 
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Figure captions 

Fig. 1. 

Half-life of ^^Cu as a function of the final spin projection in the intrinsic system Kf = K^. 
The solid lines correspond to transitions without deformation change, i. e. with f3f = f3i = 
0.35. The dashed lines correspond to = 0.35 and an oblate final shape with Pf = —0.2. 
The labels of the lines to the right of the figure are the corresponding values of Ki. 

Table. 1 

Amplitudes defined by Eq. ( ji^ j corresponding to the different Ky-values for the 9/2+ state 

at prolate and oblate shapes, respectively. 

1/2 3/2 5/2 7/2 9/2 
d% 0.901 0.405 0.089 ~- ~ 
d^P - - 0.053 0.212 0.986 
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